PREHEATING AND REHEATING IN INFLATIONARY COSMOLOGY: A 

PEDAGOGICAL SURVEY 
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Recent progress in the preheating phenomena for inflationary cosmology is reviewed. We first discuss 
estimates of the preheating time scale and particle production at the early stages of parametric 
amplification within the Mathieu and Lame approximations and we analyze their precision and 
limitations. The necessity of self-consistent calculations including the non-linearity of the field 
theory equations in an energy conserving scheme is stressed. The large A'' calculations including the 
field back-reaction are reviewed. For spontaneously broken theories the issue of symmetry restoration 
is analyzed. A discussion of the possibility and criterion for symmetry restoration is presented. (To 
appear in the Proceedings of the Paris Euronetwork Meeting 'String Gravity', Observatoire de Paris, 
June 1996). 



I. INTRODUCTION 

Research activity on inflationary cosmologies has con- 
tinued steadily since the concept of inflationary cosmol- 
ogy was first proposed in 1981 0]. 

It was recognized that in order to merge an inflationary 
scenario with standard Big Bang cosmology a mechanism 
to reheat the universe was needed. Such a mechanism 
must be present in any inflationary model to raise the 
temperature of the Universe at the end of inflation, thus 
the problem of reheating acquired further importance de- 
serving more careful investigation. The original version 
of reheating envisaged that during the last stages of infla- 
tion when the accelerated universe expansion slows down, 
the energy stored in the oscillations of the inflaton zero 
mode transforms into particles via single particle decay. 
Such particle production reheats the universe whose tem- 
perature reduced enormously due to the inflationary ex- 
pansion ||]. 

It was realized recently fl|l5|, 16 18|, that in fact, the 
elementary theory of reheating j^does not describe ac- 
curately the quantum dynamics of the flelds. 

Our programme on non-equilibrium dynamics of quan- 
tum field theory, started in 1992 is naturally poised 
to provide a framework to study these problems. The 
larger goal of the program is to study the dynamics 
of non-equilibrium processes from a fundamental field- 
theoretical description, by solving the dynamical equa- 
tions of motion of the underlying four dimensional quan- 
tum field theory for physically relevant problems: phase 
transitions out of equilibrium, particle production out 
of equilibrium, symmetry breaking and dissipative pro- 
cesses. 



The focus of our work is to describe the quantum field 
dynamics when the energy density is high. That is, a 
large number of particles per volume , where m is 
the typical mass scale in the theory. Usual S-matrix cal- 
culations apply in the opposite limit of low energy den- 
sity and since they only provide information on in — s- 
out matrix elements, are unsuitable for calculations of 
expectation values. Our methods were naturally applied 
to different physical problems like pion condensates [^,0 , 
supercooled phase transitions [^jj^ and inflationary cos- 
mology |J|,|-0. 

An analogous program has been pursued by the Los 
Alamos group, whose research focuses on non-linear ef- 
fects in scalar QED linked to the pair production in 
strong electric fields |11 , the Schwinger-Keldysh non- 
equilibrium formalism in the large N expansion [ p^ , and 
the dynamics of chiral condensates in such framework 



II. PREHEATING IN INFLATIONARY 
UNIVERSES 

As usual in inflationary cosmology, matter is described 
in an effective way by a self-coupled scalar fleld ^{x) 
called the inflaton. The spacetime geometry is a cosmo- 
logical spacetime with metric ds^ = {dt)"^ — '^(i)^ (dx)'^, 
where a{t) is the scale factor. 

The evolution equations for the fc-modes of the inflaton 
field as considered by different groups can be summarized 
as follows, 



+ M^m.)]]xk{t)^0 (2.1) 



1 



where H{t) = d{t)/a{t), and [$(.)] is the effective 
mass feh by the modes. The expression considered de- 
pends on the model (here the A model) and the ap- 
proximations made. The value of A is scenario-dependent 
but it is usually very small. [$(.)] depends on the 
scale factor and on the physical state. Therefore, it de- 
pends on the modes Xk{t) themselves in a complicated 
way. One is definitely faced with a complicated non- 
linear problem. We call 'back-reaction' the effect of the 
modes Xk{t) back on themselves through [$(.)]. 

In the initial stage, all the energy is assumed in the zero 
mode of the field (j){t) [|J|-|Jl|-|l|l . That is, the field ex- 
pectation value (/)(t) =< $(a;) >, where < • • • > stays for 
the expectation value in the translational invariant but 
non-equilibrium quantum state. For very weakly coupled 
theories, and early times, such that the back-reaction ef- 
fects of the non-equilibrium quantum fluctuations can be 
neglected, one can approximate, 
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(2.2) 



At this moment, the scale factor is set to be a constant in 



ref. |5|,y. Refs. |J|- 
from the start. In ref. [ 



consider Minkowski spacetime 
I , the scale factor is set to be a 
constant for a model without A(f>'* inflaton self-coupling. 
That is, for the classical potential llq] 



(2.3) 



one can consider as classical solution <I>(i) ~ 
$0 cos{mt), a — 0. In such case, the Mathieu equation 
approximation is exact in Minkowski spacetime. How- 
ever, the potential (2.3) is unstable under renormaliza- 
tion (a counterterm is needed from the one-loop level). 
Hence, the A = choice is a fine-tuning not protected by 
any symmetry. In a second case a massless selfcoupled 
X^"^ field in a radiation dominated universe a{t) oc \/t is 
considered in conformal time In such specific case 

the classical field equations take the Minkowski form for 

In a way or another, using th e cla ssical oscillating be- 
haviour of (/)(t), one is led by eq.(pj|) to an effective mass 
that oscillates in time. In this approximation (which 
indeed, may be very good for small coupling, see [^), 
) exhibits parametric resonance as noticed first 
p5[ |. Namely, there are allowed and forbidden 
in A;^. The modes within the forbidden bands 



eq.(2.1 



in ref 
bands 

grow exponentially whereas those in the allowed bands 
stay with bounded modulus. The growth of the modes 
in the forbidden bands translates into profuse particle 
production, the particles being created with these par- 
ticular unstable momenta. The rate of particle produc- 
tion is determined by the imaginary part of the Floquet 
index in these unstable bands. Notice that the approxi- 
mation ( |2.2[ ) breaks down as soon as many particles are 
produced. Namely, when the energy of the produced par- 
ticl es b ecomes of the order of the zero-mode energy and 
eq.(|2.2|) is no more valid. 



Now, in order to compute quantitatively the particles 
produced one needs the form of 4>{t). In ref. |p^ , p^ , p^ 
(j)(t) is approximated by a cosinus function in the calcula- 
tions. The mode equations become a Mathieu equation 
(the scale factor is set to be a constant). In ref. [T^ ] 
the Bogoliubov-Krylov approximation is used to com- 
pute estimates. In ref. p6| , p^ , estimates are obtained 
using asymptotic formulas for the Mathieu equation. In 
ref. Q the exact classical solution is used (a cn Jacobi 
function) to compute estimates. 

Let us now compare the results from the exact mode 
solutions obtained in ref. |^ with the Mathieu equation 
approximation to it. In units where = 1 and setting 
— /a 



ri{t) = \l ^ 4>{i)^ one finds 

■q{t) = 770 cn 
Vo 



k = 



(2.4) 



where cn stands for the Jacobi cosinus and we choose for 
initial conditions 77(0) = 770 , viO) = 0. 

Here r]{t) has period Auj = 4 K (k) / y^l + Tj^ , where 
K{k) is the complete elliptic integral of first kind, and 
r]{t)^ has period 2uj. 



Inserting this form for ?7(t) in eqs.(p.2|) and (2.1) yields 

Xkit) 
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(2.5) 



This is the Lame equation for a particular value of the 
coefficients that make it solvable in terms of Jacobi func- 
tions 1^^. As shown in ref. this equation has only 
one forbidden band for positive k^ going from fc^ = to 

k"^ = One can choose Floquet solutions of eq.(pT5|) 
that fullfil the relation 



C/fc(i + 2c^) = e'^W Uk{t), 



(2.6) 



where the Floquet indices F{k) are independent of t. In 
the forbidden band the F(k) posses an imaginary part. 
The production rate is determined by the imaginary part 
of the Floquet index. 

The exact form of F(k) results ||], 

F{k) = -2iK{k) Z{2K{k) v)+t: 

where Z{u) is the Jacobi zeta function p5[ and w is a 
function of k in the forbidden band defined by 

k^ ^CTL{2K{k)v,k) , Q <v <]^. (2.7) 

All these elliptic functions posses fastly convergent ex- 
pansions in powers of the elliptic nome 

q ^ ^-rrK'{k)/KCk) _ 
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Since < fc < l/\/2 [see eq.(2^)], we have 

< g < ^ 0.0432139. .. . (2.8) 

Then, 

F{k) = 4j TT g sin(27rw) [l + 2q cos 2ttv + 0(g^)] + tt . 

(2.9) 

The imaginary part of this function has a maximum at 
fc = fci = i ?7o (1 - g) + 0(g2) where g 



T = ImF(ki) = 4 7rg + 0(g^) 



(2.10) 



This simple formula gives the maximum of the imaginary 
part of the Floquet index in the forbidden band with a 
precision better than 8. 10~^. q can be expressed in terms 
of ?7o as follows H 



2 (1 + 772)1/4 + (l + ,y2/2)i/4 



with an error smaller than ~ 10^^. 

Let us now proceed to the Mathieu equation analysis 
of this problem. The cn Jacobi function can be expanded 
as pi 



width of the band is well approximated but not its abso- 
lute position. The numerical values of the maximum of 
the imaginary part of the Floquet index are given in Ta - 
ble I and compared with the exact values from eq.( 2.10 ). 

We see that the Mathieu approximation underesti- 
mates the exact result by a fraction ranging from 13% 
to 39%. The second forbidden band in the Mathieu 
equation yields ^Mathieu = 0.086 . . . for 770 — > 00. This 
must be compared with Thame — (no further forbidden 
bands). 

In ref. p7[ |, an even larger discrepancy between Lame 
and Mathieu Floquet indices has been reported within a 
different approximation scheme. 

It must be noticed that the Floquet indices entering in 
exponents produces a dangerous error propagation. For 
example, the number of particles produced during re- 
heating is of the order of the exponential of 2T times the 
reheating time in units of tt/ 1 -I- 77Q . An error of 25% in 
T means an error of 25% in the exponent. (For instance, 
one would find 10^ instead of 10^^). 

The mode equations ( |2.l[ ) apply to the self-coupled 
A scalar field. Models for reheating usually contain 
at least two fields: the inflaton and a lighter field <j(x) 
in which the inflaton decays. For a gcP'^^ coupling, the 
mode equations for the a field take the form 



6 18| 



cn(z, fc) = (1 — q) cos(l — 4g)z -f g cos 32; -|- O(q^) 
To zeroth order in q we have 



1 + cos{2t^l + ril) 



0{q) 



and 2lo = IT / \/\ + riQ+0{q). Under such approximations 
eq.(2.5) becomes the Mathieu equation 



+ (a — 2(7Cos2z) y{z) = 



(2.11) 



Vk+3H{t)Vk- 



|F[$(.),a(.)])y,.(t) 







(2.12) 



A new dimensionless parameter ^ appears here. Neglect- 
ing the a and $ backreaction, we have 



(2.13) 



In ref. [|6|J^ , it is shown that abundant particle produc- 
tion (appropriate for reheating) shows up even for g = \. 



where 



a = 1 



fc2 _ !k 
2^ 2_ 

77o' + l 



4(772 + 1) 



and z = yj-ql + 1 t. Notice that < g < 1/4 in the 
present case. Eq.(2.11) posses an infinite number of for- 
bidden bands for W> 0. The lower and upper edge of 
the first band are respectively p^] 



fc2 



1 - 



4 



tiU 



25(772 + 1) 



210(772 
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2i"(77o2 



1) 



770 


T Lame 


Math-ieu 


%error 


1 


0.2258 . . . 


0.20... 


13% 


4 


0.4985 . . . 


0.37... 


35% 


770 — » CO 


^■Ke"" = 0.5430 . . . 


0.39... 


39% 



TABLE I. The maximum of the imaginary part of the 
Floquet index T for the Lame equation and for its Mathieu 
approximation. 



These values must be compared with the exact result 
for the Lame equation (|2.5|) : fc? f = 



fc2 

sup 



f . The 
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Eqs.( ^.12 )-( 2.13| ) becomes a Lame equation when ri(t) 
is approximated by the classical solution in Minkowski 
spacetime given by {^A\). Such Lame equation is solvable 
in closed form when the couplings g and A are related as 
follows M 



2l 
A 



n{n + 1) , n = 1, 2, 3, , 



In those cases there are n forbidden bands for fc^ > 0. 
The Lame equation exhibits an infinite number of for- 
bidden bands for generic values of j. The Mathieu 
and WKB approximations has been applied in the non- 
exactly solvable cases [p^|~p^. However, as the above 
analysis shows (see Table I) such estimations cannot be 
trusted quantitatively. The only available precise method 
consists on accurate numerical calculations as those of 
ref. H"!^ (where the precision is at least 10"^). 

Estimates in the cosinus approximation for FRW- 
de Sitter backgrounds and open universes using the 
Bogoliubov-Krylov approximation are given in ref. 
In ref. [ p2[ the Bogoliubov-Krylov approximation is ap- 
plied to the large N equations. 

Applications of preheating to various relevant aspects 
of the early cosmology are considered in ref. |23]. 



As soon as the quantum fluctuations grow and cease 
t o be negligible compared with the the classical piece 
(2.13), all the approximations discussed so far (Lame, 
Mathieu, etc.) break down. This time is the so-called 
preheating time trek H ■ O ne can estimate trek by equat- 



ing the zero mode energy ( p. 131) with the estimation of 
the quantum fluctuations derived from the unstable Flo- 
quet modes [|l^,P . Such estimation yields when the Lame 
Floquet indices are used M, 



t. 



reh 



1 iV(l 



%V2) 



(2.14) 



where 



B = 8./l + 77o'9(l-4g) + 0(g^) , 



N = 



(4 + 3r;g)74T5^ 
r7§(l + ry2)3/4 



[l + 0{q)] . (2.15) 



Where B is determined by the maximum of the imagi- 
nary part of the Floquet index in the band. It is now clear 
that a few percent correction to the Floquet indices will 
result in a large error in the estimate of the preheating 
time scale. 



III. THE END OF PREHEATING AND BEYOND: 
SELF-CONSISTENT LARGE TV CALCULATIONS 

In order to compute physical magnitudes beyond trek, 
one must solve self-consistently the field equations in- 
cluding the back reaction. In ref. this is done for 
the N ^ oo limit and in ref. m to one-loop order. 



In the large N limit the zero mode and fc-modes renor- 
malized equations take the form, 



ii + v + v^' + g v{t) s(t) - , 



Vk{t) = 



(3.1) 
(3.2) 



where (/S(i) is given by 



9{k-K) 

2fc3~ 



[-fll + ri\t)+gm] 



(3.3) 



We choose the initial state such that a,t t — the quan- 
tum fluctuations are in the ground state of the oscillators. 
That is. 



1 



r){0) = 770 , viO) = . 

In the one-loop approximation the term g E(t) is absent 
in eq.(n 

Eqs (|3.2D-(^) were generahzed to cosmological space- 
times (including the renormalization aspects) in ref. [Q. 

Eqs.( |3.2[ ) is an infinite set of coupled non- linear dif- 
ferential equations in the fc-modes ipk{t) and in the zero 
mode r]{t). We have numerically solved egs. _(|3.2[ ) for a 
variety of couplings and initial conditions ||6|-|8| . In figs. 
1-3 we display 77 (i), S(t) and the total number of created 
particles N{t) as a function of t for 770 = 4 and g = 10~^^ 
One sees that the quantum fluctuations are indeed 
very small till we approach the rehea ting time treh- As 
one can see from the figures, eq.( 2.14 ) gives a very good 
approximation for treh and for the behaviour of the quan- 
tum fluctuations for t < treh- For times earlier than treh, 
the Lame estimate for the envelope of S(t) and N{t) take 
the form §, 



Ne. 



At) 



„Bt 



N yft 



^est — env it) . (3.4) 



Stt^ ^4 + 5772 



where B and N are given by eq. ( 2.15 ) . 

This analysis presents a very clear physical picture of 
the properties of the "gas" of particles created during the 
preheating stage. This "gas" turns to be out of thermo- 
dynamic equilibrium as shown by the energy spectrum 
obtained, but it still fulfills a radiation-like equation of 
state p ~ e/3 ||]. 

In approximations that do not include the backreaction 
effects there is infinite particle production, since these 
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do not maintain energy conservation The full 

backreaction problem and the approximation used in our 
work exactly conserves the energy. As a result, parti- 
cle production shuts-ofF when the back-reaction becomes 
important ending the preheating stage. 

In ref. [Q results of ref. [D are rederived with a differ- 
ent rcnormalization scheme. 



IV. BROKEN SYMMETRY AND ITS QUANTUM 
RESTORATION 

Up to now we discussed the unbroken symmetry case 
= 1 > 0. In the spontaneously broken symmetry case 
771^ = — 1 < 0, the classical potential takes the form 



If 



(4.1) 



In refs. ||6|,g| the field evolution is numerically solved in 
the large N limit for ?7t,^ < 0. Analytic estimations are 
given in ref. Q for early times. 

In ref. |^ it is shown that for small couplings and small 
r]Q, the order parameter damps very quickly to a constant 
value close to the origin. That is, the order parame- 
ter ends close to the classical false vacuum and far from 
the classical minimum r/ = 1. The fast damping is ex- 
plained, in this spontaneously broken symmetry case by 
the presence of Goldstone bosons. The massless Gold- 
stone bosons dissipate very efficiently the energy from 
the zero mode. We plot in figs. 4 and 5, ri{t) and 
for 7/0 = 10^^ and g = 10~^^ |8|. The profuse particle 
production here is due to spinodal unstabilities. 

In ref. [^, it is claimed that quantum fluctuations 
have restored the symmetry in the situations considered 
in ref. However, the presence of massless Goldstone 
bosons and the non-zero limiting value of the order pa- 
rameter clearly show that the symmetry is always broken 
in the cases considered in ref. |^,|| . These results rule out 
the claim in ref. |l9| ]. In particular, the asymptotic value 
of the order parameter results from a very detailed dy- 
namical evolution that conserves energy. 

In the situation of 'chaotic initial conditions' but with 
a broken symmetry tree level potential, the issue of sym- 
metry breaking is more subtle. In this case the zero mode 
is initially displaced with a large amplitude and very 
high in the potential hill. The total energy density is 
non-perturbatively large. Classically the zero mode will 
undergo oscillatory behavior between the two classical 
turning points, of very large amplitude and the dynam- 
ics will probe both broken symmetry states. Even at the 
classical level the symmetry is respected by the dynam- 
ics in the sense that the time evolution of the zero mode 
samples equally both vacua. This is not the situation 
that is envisaged in usual symmetry breaking scenarios. 
For broken symmetry situations there are no finite energy 
field configurations that can sample both vacua. In the 
case under consideration with the zero mode of the scalar 



field with very large amplitude and with an energy den- 
sity much larger than the top of the potential hill, there is 
enough energy in the system to sample both vacua. (The 
energy is proportional to the spatial volume). Paramet- 
ric amplification transfers energy from the zero mode to 
the quantum fluctuations. Even when only a fraction of 
the energy of the zero mode is transferred thus creating 
a non-perturbatively large number of particles, the en- 
ergy in the fluctuations is very large, and the equal time 
two-point correlation function is non-perturbatively large 
and the field fluctuations arc large enough to sample both 
vacua. The evolution of the zero mode is damped because 
of this transfer of energy, but in most generic situations 
it does not reach an asymptotic time-independent value, 
but oscillates around zero, sampling the tree level min- 
ima with equal probability. This situation is reminiscent 
of finite temperature in which case the energy density 
is finite and above a critical temperature the ensemble 
averages sample both tree level vacua with equal prob- 
ability thus restoring the symmetry. In the dynamical 
case, the "symmetry restoration" is just a consequence 
of the fact that there is a very large energy density in 
the initial state, much larger than the top of the tree 
level potential, thus under the dynamical evolution the 
system samples both vacua equally. This statement is 
simply the dynamical equivalent of the equilibrium finite 
temperature statement that the energy in the quantum 
fluctuations is large enough that the fluctuations can ac- 
tually sample both vacua with equal probability. 

Thus the criterion for symmetry restoration when the 
tree level potential allows for broken symmetry states is 
that the energy density in the initial state be larger than 
the top of the tree level potential. That is when the am- 
plitude of the zero mode is such that ^(770) > V{0). In 
this case the dynamics will be very similar to the unbro- 
ken symmetry case, the amplitude of the zero mode will 
damp out, transferring energy to the quantum fluctu- 
ations via parametric ampliflcation, but asymptotically 
oscillating around zero with a fairly large amplitude. 

To illustrate this point clearly, we plot in flgs. 6 
and 7, ri{t) and I](t) for rjo = 1.6 > \/2 [and hence 
V{t]o) > V{0), see eq.{^)] and g = IQ-^. We find the 
typical behaviour of unbroken symmetry. Notice again 
that the effective or tree level potential is an irrelevant 
quantity for the dynamics, the asymptotic amplitude of 
oscillation of the zero mode is rj « 0.5, which is smaller 
than the minimum of the tree level potential rj = 1 but 
the oscillations are symmetric around rj = 0. 

Since the dynamical evolution sampled both vacua 
symmetrically from the beginning, there never was a 
symmetry breaking in the first place, and "symmetry 
restoration" is just the statement that the initial state 
has enough energy such that the dynamics probes both 
vacua symmetrically despite the fact that the tree level 
potential allows for br oken symmetry g round states. 

In their comment ( hcp-ph/960834l| ) , KLS seem to 
agree with our conclusion that in the situations stud- 
ied in our articles pM|, when the expectation value is 
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released below the potential hill, symmetry restoration 
does not occur. 

We will present a deeper analytical and numerical 
study of the subtle aspects of the dynamics in this case 
in a forthcoming article M. 



V. LINEAR VS. NONLINEAR DISSIPATION 
(THROUGH PARTICLE CREATION) 

As already stressed, the field theory dynamics is un- 
avoidable nonlinear for processes like preheating and re- 
heating. It is however interesting to study such processes 
in the amplitude expansion. This is done in detail in refs. 
[^,0. To dominant order, the amplitude expansion means 
to linearize the zero mode evolution equations. This ap- 
proach permits an analytic resolution of the evolution 
in closed form by Laplace transform. Explicit integral 
representations for ri{t) follow as functions of the initial 
data 1^,01 . Moreover, the results can be clearly described 
in terms of S-matrix concepts (particle poles, production 
thresholds, resonances, etc.). 

Let us consider the simplest model where the inflaton 
$ couples to a another scalar a and to a fermion field ip, 
and potential Q 

In the unbroken symmetry case (rn% > 0) the inflaton 
is always stable and we found for the order parameter 
(expectation value of $) evolution in the amplitude ex- 
pansion 1^], 



v{t) = 



-1 OS(im^) 

-L — a ^5 

dm- 



cos[m$i] 



wl]/(ijj) COS ujt duj 



(5.1) 



where Ephysicai(jw±0^) = S]fl(w) ±iS]/(cij) is the inflaton 
self-energy in the physical sheet, rji = 77(0) and //(O) = 
0. The first term is the contribution of the one-particle 
pole (at the physical inflaton mass). This terms oscillates 
forever with constant amplitude. The second term is the 
cut contribution r]{t)cut corresponding to $ ^ $ -|- 2a. 
In general, when 



S;/(a; — > ^threshold) 



^^^thrcshold 



the the cut contribution behaves for late times as 
2rii B Wthroshoid r(l -I- a) 



S_R,(wthr 



sholdJ 



l^thrcsholdi + U 



(5.2) 



Here, Wthioshoid 
loops, § 



S/(a>) 



771$ + 2Ma- and a = 2 since to two- 



(47r)4 (to$ + 2M^)7/2 



[uj^-{m^+2Mcr) 



F = 



1 



A=l 



2i2 



In the broken symmetry case (m% < 0) we may have 
either M < 2ma or M > where M is the physical 
inflaton mass. ( M = |m$|\/2 at the tree level). In the 
first case the inflaton is stable and eq.(^|^) holds. How- 
ever, the self-energy starts now at one-loop and vanishes 
at threshold with a power a = 1/2. For M > m„ the 
inflaton dccomes a resonance (an unstable particle) with 
width (inverse lifetime) 



8ttM V ' M2 
This pole dominates rjit) for non asymptotic times 

S{t) c^S.A e~^*/2 cos{Mt + 7) , (5.3) 

where 



7 = -^^ (54) 



In summary, eq.(5.3) holds provided: a) the inflaton is 
a resonance and h) t < F^^ln(F/MCT). For later times 
the fall off is with a power law t~'^/^ determined by the 
spectral density at threshold as before 0. 

In ref. the selfconsistent nonlinea r ev olution is com- 
puted to one-loop level for the model (5.1). In fig. 8 r]{t) 



is plotted as a function of time for A = 5 = I.Gtt , y — 
0, = 0.2771$, 77(0) = 1 and 77(0) = 0. 

Figure 8 shows a very rapid, non-exponential damp- 
ing within few oscillations of the expectation value and 
a saturation effect when the amplitude of the oscillation 
is rather small (about 0.1 in this case), the amplitude re- 
mains almost constant at the latest times tested. Figures 
8 and 9 clearly show that the time scale for dissipation 
(from fig. 8) is that for which the particle production 
mechanism is more efficient (fig. 9). Notice that the to- 
tal number of particles produced rises on the same time 
scale as that of damping in fig. 8 and eventually when 
the expectation value oscillates with (almost) constant 
amplitude the average number of particles produced re- 
mains constant. This behaviour is a close analog to the 
selfcoupled inflaton for unbroken symmetry (fig.l). The 
amplitude expansion predictions are in qualitative agree- 
ment with both results. 

These figures clearly show that damping is a conse- 
quence of particle production. At times larger than about 
40 777 J ^ (for the initial values and couplings chosen) there 
is no appreciable damping. The amplitude is rather small 
and particle production has practically shut off. If we had 
used the classical evolution of the expectation value in the 
mode equations, particle production would not shut off 
(parametric resonant amplification), and thus we clearly 



6 



see the dramatic effects of tfie inclusion of tfie back reac- 
tion. 

In ref. Q the broken symmetry case to|, < is then 
studied. Figures 11-13 show 77(t) vs t, Ma{T) vs t and 
Nq^a{T = 200) VS q respectively, for A/Stt^ = 0.2; g/\ = 
0.05; m„ = 0.2|m<t,|; 77(0) = 0.6; 77(0) = 0. 
Notice that the mass for the linearized perturbations 
of the $ field at the broken symmetry ground state 
is \/2|w$| > 2m„. Therefore, for the values used in 
the numerical analysis, the two-particle decay channel 
is open. For these values of the parameters, linear re- 
laxation predicts exponential decay with a time scale 
Tr-ei ~ 300 (in the units used). Figure 11 shows very 
rapid non-exponential damping on time scales about six 
times shorter than that predicted by linear relaxation. 
The expectation value reaches very rapidly a small am- 
plitude regime, once this happens its amplitude relaxes 
very slowly. In the non-linear regime relaxation is clearly 
not exponential but extremely fast. The amplitude at 
long times seems to relax to the expected value, shifted 
slightly from the minimum of the tree level potential at 
Tj = 1. This is as expected from the fact that there are 
quantum corrections. Figure 12 shows that particle pro- 
duction occurs during the time scale for which dissipa- 
tion is most effective, giving direct proof that dissipation 
is a consequence of particle production. Asymptotically, 
when the amplitude of the expectation value is small, 
particle production shuts off. We point out again that 
this is a consequence of the back- reaction in the evolution 
equations. Without this back-reaction, as argued above, 
particle production would continue without indefinitely. 
Figure 13 shows that the distribution of produced par- 
ticles is very far from thermal and concentrated at low 
momentum modes k < |m$|. This distribution is quali- 
tatively similar to that in the unbroken symmetry case, 
and points out that the excited state obtained asymptot- 
ically is far from thermal. 

In ref. |Q the case where the inflaton is only coupled 
to fermions is studied {g — 0, y ^ 0). The damping 
of the zero mode is very inefficient in such case due to 
Pauli blocking. Namely, the Pauli exclusion principle 
forbids the creation of more than 2 fermions per momen- 
tum state. Pauli blocking shuts off particle production 
and dissipation very early on. 

VI. FUTURE PERSPECTIVES 

The preheating and reheating theory in inflationary 
cosmology is currently a very active area of research 
in fast development, with the potential for dramatically 
modifying the picture of the late stages of inflationary 
phase transitions. 

As remarked before, estimates and field theory cal- 
culations have been done mostly assuming Minkowski 
spacetime. Results in de Sitter [|l^ and FRW back- 
grounds pql are just beginning to emerge. A further 



important step will be to consider the background dy- 
namics. Namely, the coupled gravitational and matter 
field dynamics. The matter state equations obtained in 
Minkowski |^ and de Sitter backgrounds ^ give an indi- 
cation, through the Einstein-Friedmann equation, of the 
scale factor behaviour. 
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Figure 11: The inflaton in the broken symmetry case 
coupled to a hghter scalar a. rj{T) vs r for the values 
of the parameters y = 0; A/Stt^ = 0.2; g = X; = 
0.2|m^|; ?7(0)=0.6; ?7(0) = 0. 
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Figure 12: A/'ct(t) vs. t for the same value of the 
parameters as fig. 11. 
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Figure 13: Mq^ai^ = 200) vs. q for the same values 
as in fig. 11. 
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